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The functional renormalisation group is applied to the effective action for scatter- 
ing of two nonrelativistic fermions. The resulting physical effective action is shown to 
contain the correct threshold singularity. The corresponding "bare" action respects 
Galilean invariance up to second order in momenta. Beyond that order it contains 
terms that violate this symmetry and, for the particular regulator considered, non- 
analytic third-order terms. The corresponding potential can be expanded around a 
nontrivial fixed point using the power counting appropriate to a system with large 
^ ' scattering length. 

■ The functional or "exact" renormalisation group (RG) is proving to be a powerful tool 

\ for the study of phase transitions in a variety of systems, from condensed matter to particle 



■ physics [l|, l2|. Recently there has been particular interest in applying it to pairing in dense 



Q 

matter consisting of nonrelativistic fermions, such as nuclear matter or cold atoms in traps 
^ ! Q, S]- In this context, it has been used to examine the cross-over between BCS superfluidity 



O 



and Bose-Einstein condensation of tightly bound pairs. 



^ \ In this approach, a scale-dependent regulator is used to interpolate between a bare action 

' and the physical effective action. This regulator is designed to suppress the contributions of 
fiuctations with momenta smaller than the cut-off scale, k. As k ^ 0, it vanishes so that all 
fluctuations are integrated out and we arrive at the full effective action describing physical 
scattering. The starting point for the RG flow is a suitably parametrised "bare" action at 
some large starting scale K. This is often expressed as an expansion in derivatives of the 
fields. To fix the parameters in it, we can evolve down to the physical limit in vacuum and 
there we can match, for example, to the two-body scattering matrix. 

The input into these calculations of fermionic matter is thus in the form of two-body (and 
ultimately also three- and more-body) forces between the particles in vacuum. This means 
that there is a close connection to the effective field theories that are being developed for 



few-body systems, especially in nuclear physics and, in particular, to the Wilsonian 

■iG methods that have been developed to determine the power countings in these theories 
8]. In fact, the Kyushu group has applied the functional RG to nucleon-nucleon scattering 
91] and shown that the results agree with the Wilsonian approach. A first application to 
three-body systems has been made more recently jlO |. 

Here I apply the "exact" RG to the scattering of two nonrelativistic fermions, using the 
version based on the Legendre-transformed effective action and introducing a boson 

field to describe correlated pairs of fermions ("dimers") as in Refs. js, 0, sl- In the present 
work I do not make a derivative expansion of the boson self-energy but instead I keep terms 
to all orders in energy and momenta. This provides an extension of the calculations in 
Ref. ^ and Sec. IV of [lO] that allows both an analysis of the scaling properties, as in 
Ref. 81], and the identification of artefacts caused by the regulator, such as violations of 
Galilean invariance. 

In describing systems of fermions with attractive forces, it is convenient to represent the 
interactions in terms of boson fields. These fields describe correlated, but not necess arily 



bound, pairs of fermions, like the dibaryon fields often used in nuclear effective theories |ll ]. 
Here I consider a single type of fermion, with two spin states, and I introduce one boson 
field to describe the interactions between them. The Fourier transforms of these fields are 
denoted by ip{q) and (j){q), respectively. The starting point is the following ansatz for the 
effective action F, 

,2 



F[7A,^^0,0^;fc] = / d\ 



2M 

~^(2^ / ^^^^^^^2 [i^l){.qi + q2Y^{,q2f(72ij{,qi) 

V^(gi) V2V^(g2)^^0(gi + g2)) + ■ ■ ■ • (1) 

This is essentially the same as the action considered in Refs. , except for the inverse 

boson propagator H(go, I'-, k) which I do not expand in powers of energy or momentum. 

In Eq. ([T]), I have chosen to express the interactions between the fermions entirely in terms 
of the boson field. The action thus does not include any explicit two-body interactions and 
one can regard as an auxiliary field arising from bosonisation of these terms. If I make 
this choice at the starting scale K then, at least in vacuum, it is maintained during the 
RG evolution; there are no one-particle-irreducible loop diagrams that generate two-body 
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interactions between the fermions. In vacuum, there are also other simphfications. In 
particular, conservation of fermion number means that there are no boson loops that dress 
either the fermion propagator or the fermion-boson coupling and so these retain their free 
forms. 

The action evolves in the usual way as the regulator scale is changed [l, 2|: 



dkT = +- Tr 



FF. 



Tr 



B 



T(2) _ R)-l 



BB 



(2) 



where F*^^-* is the matrix of second derivatives of the action with respect to the fields. In 
vacuum we shall need only the fermionic part, which involves 



(r(2) - r) 



FF 



2M 



RFiq,k) ig(j)a2 



\ 



-ig(f)^a2 qo + — + RF{q,k) 



(3) 



Note that, although no approximation has yet been made, the driving terms on the right 
hand side of Eq. ([2]) have the form of one-loop integrals. 
The evolution of the boson self-energy can be found from 

52 



dkIl(Po, P, k) 
where, in vacuum, only fermion loops contribute: 

<t>=0 



(4) 



50(P) 5(p{P) 



+i Tr 



<t>=0 



lT-(3)t /-n(2) 



^ r^]^^ - Rf) 



(5) 



Here the necessary third derivatives of F are 



.(3) 

FF(I) 



d 



.(2) 
FF 




(6) 



and the matrix trace gives 



tr 



^(2) 



RF)-^Fl^LfF?) 



{SkKp) \^ FF ~ ^^F) FF4}\^ FF ~ ^'-F) ^FF<pK-^FF 



Rf)-^f?) 



^(2) 



Rf 



dkRF{q, k) 



1 2 



qo + EpRiq, k) - ie go + ^0 - Epniq + P,k)+ ie 



(7) 
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where 

EpR^q. k) = q^/2M + Rp^q, k). (8) 

After performing the contour integral over go and symmetrically routing the external 
momentum P (by shifting q q — P /2 and using the symmetry under P —P)i the 
driving term can be written in the form of a derivative with respect to k: 

d?q 1 



5(p{P) 5(p{Py 



= g dk 

J7=0 



(27r)3 EpRiq - P/2, k) + EpR^q + P/2, k) - Po - ie 

(9) 

This allows us to relate n(Po5 P] k) to the physical self-energy (at A; = 0), by integrating the 
evolution equation to get 



n(Po,P;A;) - n(Po,P;0) 



d?q r 1 



(27r)3 [ EFR{q - P/2, 0) + EpRiq + P/2, 0) - Pq - ie 

1 

~ EFR{q - P/2, k) + EpRiq + P/2, k) - Po - ie J ' 

(10) 

(From now on, all expressions are three-dimensional and P denotes the magnitude of P.) 
A convenient choice of regulator is the one used by Diehl et al. jsl, 

RF{q,k) = ^^e{k-q). (11) 

This is just the nonrelativistic version of the "optimised" cut-off for boson fields introduced 
by Litim [l2|. It is a variant of the sharp cut-off which has the advantage of leading to 
a momentum-independent fermion propagator for three-momenta smaller than k. With it, 
the integrand in Eq. ( fTOl) vanishes outside the intersecting spheres, \q ± P/2| < k. In the 
region between the spheres, the integrand for g^ > is 

M 2M 
g2 + g2 + P74 - MPo ~ g^ + (g, + P/2)2 + P - 2MPo ' ^^^^ 

where I have chosen the 2;-axis along the direction of P and I have defined g^ = ^Jq'^ + g^. 
It is convenient to integrate this over the union of the two spheres, < g^ < A; + P/2, 
< g^ < Vfc2-(g^-P/2)2. The result can be corrected by adding the integral of 

2M M 



ql + (g, + P/2)2 + fc2 - 2MPo k'^ - MPq ' 



(13) 
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over the intersection of the spheres, < < A; — P/2, < < ^Jk'^ — {qz + P/2y. The 
symmetry of the integrand means that the regions with < make identical contributions 
to the corresponding ones with gz > 0. 

After integrating over q, the resulting form of the self-energy is 

u{Po,P;k) - n(Po,P;0) 



- v^MPo -i^V4 In 



'k + P/2 + v^MPo -PV4' 
k + P/2- JMPo-P^/4: 



1 



A;2 - MPn 



10 3 P^ 

— A;^ - 4A;MPn --k^P + 2MPnP 



+ 4y/k^ -2MPo 



arctan 



2 

fc + P 



24 



— arctan 



VA;^ - 2MPo 
k 



P2 - MP 



In 



_ 2MPo 

k'^ + kP + P^/2 - MPo 
A;2 - MPo 



(14) 



Here I have assumed that the scale k is large {k > y/2MP^). 

The on-shell fermion-fermion scattering amplitude is directly related to the boson self 
energy in the physical limit {k 0) by 

1 1 



n(Po,P; 0), 



(15) 



Tip) 

where p = \/ M Pn — P^ /4 is the relative momentum of the two fermions. The effective-range 
expansion [131] can be used to write this amplitude as 



T{p) 



M 

— \-ip-- + -r^p^ + 
An \ a 2 



(16) 



Comparing this with the result f[T^ . we see that the RG evolution generates the correct 
nonanalytic term describing the threshold (oc i a/ MPq — P2/4) . As a result this cancels to 
leave a self-energy n(Po, P; k) that is real and can be expanded as a power series in energy 
and momentum for large k. 

Inserting the physical amplitude into Eq. (HM and expanding in powers of Pq and P, we 
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find that Il{Po, P; k) can be written 



U{Po,P; k) 



47r2 



4 , TT 

6 a 



« TT 

3fc ~ 2 



(17) 



The use of a regulator that does not respect Galilean invariance means that the evolution 
generates terms that are not invariant. To cancel these and leave the correct physical 
amplitude at = 0, the self-energy at the starting scale k = K must also contain such 
terms. The leading term that violates Galilean invariance is of order P^, and others appear 
at all higher orders. 

The order-P^ term is also nonanalytic in P^ and so could not be generated by a power 
of V^. Such terms are artefacts of regulators like the sharp cut-off used here, Eq. f|TT]) . As 
discussed in Refs. jl4, 



151, 



16l . they arise from the nonlocalities introduced by the disconti- 

nn 

nuity at the surface q = k. The cut-offs used here and in Refs. [5|, 112|] vanish at this surface 
and so the first such term appears at order P^ A cut-off with a stonger discontinuity, 
for example 



RF{q,k) 



2M 



0{k-q), 



(18) 



would lead to nonanalytic dependence appearing at order P 

As in Ref. jsl, we can define a boson wave-function renormalisation factor by 

d 



Z^k) 

and a kinetic-mass renormalisation by 

Zm{k) = 



dPn 



n(Po,P; k) 



Po=P=0 



d 



4M 



Qp2 



n(Po,P;A;) 



Po=P=0 

From the expansion of n(Po, P; k) above, we find that these are 



Z^{k) = Zm{k) 



g^M^ f 8 



TT 



(19) 



(20) 



(21) 



47r2 \?,k 2 

The fact that Z^{k) and Zm{k) are identical is a consquence of the Galilean invariance of 
the vacuum (which is respectec. at _this order). This will not be the case for calculations of 



fermionic matter as in Refs. [3|, 



3, [s], where there is a preferred rest frame. These factors 
diverge as /c as a result of the nonanalytic energy-dependence of the threshold term 
in the physical amplitude. This provides a warning that a derivative expansion may not 
be valid for the physical effective action, even it can be used to parametrise the input bare 
action. 
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Evaluated at the starting scale the expressions (12T!) can be used to provide initial 
conditions on the RG evolution. They show that the physical effective range is encoded in 
the initial values for these renormalisation factors. Note that their initial values can become 
negative if the effective range is positive and we try to start from a scale K > IG/Svrre 
(i.e. above the scale of the underlying physics). 

For a starting scale K larger than the momenta of interest, we can relate Il{Po, P', K) 



to an effective potential as in Ref. jo]. Allowing for violations of Galilean invariance, this 
potential is defined by 

= \u{{p' + Py4)/M,P-K) 

M|4_7r/8 ttAo 1^3 



V{p,P-K) g 



,^--K+-+\ rAp' pij^...\, (22) 

47r2 1 3 a \m 2 24:K^ J ^ ^ 

This can then be used, following Ref. 8|, to define a rescaled potential V = {MK /27c'^)V 
and dimensionless momentum variables p = p/K, and P = P/K. From this we see that the 
part of V coming from the integral in Eq. f|T^ is independent of K. As already noted in 
the conext of a derivative expansion [9] of the action, this is the nontrivial fixed point of the 
RG, corresponding to a system with infinite scattering length. The different regulator means 
that detailed form of the fixed point differs from that found with a simple sharp cut-off , 
but the fact that it contains a constant term of order unity (in this case —2/3) is universal. 
Similarly, the terms in the potential from the effective-range expansion scale with the same 
powers of K as found in Ref. 8|. The resulting power counting is the usual one for systems 
with large scattering lengths 6|. 

This calculation of the untruncated two-body effective action illustrates several features 
of more general significance for applications of the functional RG to few-body systems and 
fermionic matter. The lack of Galilean invariance of the regulator means that non-invariant 
terms must be included in the initial bare action. With a sharp cut-off like the one used 
here, the bare actions will also need to contain terms with nonanalytic dependences on 
P^, which would not be present in a standard derivative expansion. The need for such 
unphysical terms, with carefully tuned coefficents, is a potentially significant complication 
for applications of this approach to matter. It has not had to be addressed so far, but only 
because the actions used have been truncated at too low orders for these effects to arise 
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The RG evolution generates the correct square-root singularity associated with the thresh- 
old for two-body scattering. However, although this is easy to identify in the analytic solu- 
tion of the RG equations presented here, the corresponding singularities will require careful 
treatment as the physical limit is approached in numerical integrations. For example, if the 
action is expanded in powers of derivatives, the resulting boson wave-function renormalisa- 
tion diverges as A; ^ 0. 

Finally, the two-body action for large scales can be expressed in terms of an effective 
potential. This can be expanded around a nontrivial fixed point using the power counting 
appropriate to a system with large scattering length. The terms in the boson self-energy are 
then in one-to-one correspondence with the terms in the effective-range expansion. 
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